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ABSTRACT 

From observations collected with the ESPaDOnS spectropolarimeter, we report the 
discovery of magnetic fields at the surface of the mildly accreting classical T Tauri 
star (cTTS) V2129 Oph. Zeeman signatures are detected, both in photospheric lines 
^ . and in the emission lines formed at the base of the accretion funnels linking the disc to 

QLh- the protostar, and monitored over the whole rotation cycle of V2129 Oph. We observe 

that rotational modulation dominates the temporal variations of both unpolarised and 
^ . circularly polarised line profiles. 

^ ■ We reconstruct the large-scale magnetic topology at the surface of V2129 Oph 

^ \ from both sets of Zeeman signatures simultaneously. We find it to be rather complex, 

with a dominant octupolar component and a weak dipole of strengths 1.2 and 0.35 kG 
respectively, both slightly tilted with respect to the rotation axis. The large-scale field 
is anchored in a pair of 2 kG unipolar radial field spots located at high latitudes and 

■ coinciding with cool dark polar spots at photospheric level. This large-scale field ge- 
ometry is unusually complex compared to those of non-accreting cool active subgiants 
with moderate rotation rates. 

As an illustration, we provide a first attempt at modelling the magnetospheric 
topology and accretion funnels of V2129 Oph using field extrapolation. We find that 
the magnetosphere of V2129 Oph must extend to about 7 to ensure that the 

■ footpoints of accretion funnels coincide with the high-latitude accretion spots on the 
I stellar surface. It suggests that the stellar magnetic field succeeds in coupling to the 

accretion disc as far out as the corotation radius, and could possibly explain the 
slow rotation of V2129 Oph. The magnetospheric geometry we derive qualitatively 
. reproduces the modulation of Balmer lines and produces X-ray coronal fiuxes typical 

of those observed in cTTSs. 
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1 INTRODUCTION 

T Tauri stars (TTSs) are young low-mass stars that have 
emerged from their natal molecular cloud core. Among 
them, classical TTSs (cTTSs) are those still surrounded by 
accretion discs. CTTSs are strongly magnetic; their fields 
are thought to be responsible for disrupting the central re- 
gions of their accretion discs and for channelling the disc 
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material towards the stellar surface along discrete accre- 
tion funnels. This magnetically-channelled accretion can 
deter mine th e angul ar momentum e volution of protostars 



(e.g 



Konid 1991; S hu et al.1 [Tiii : ICameron Campbelll 
Romanova et al. 2004) as well as their internal struc- 
ture. Studying magnetospheric accretion processes of cTTSs 
through both observations and simulations is thus crucial for 
our understanding of stellar formation. 

The Zeeman broadening of unpolarised sp ectral lines 
, IJohns-Kru h et al.1 Il999: 'Ya ng et al.1 12005). has been 
very successful at determining mean surface magnetic field 
strengths on cTTSs, even for complex field topologies, but 
is rather insensitive to the large-scale structure of the 
field. Spectropolarimetric techniques, aimed at detecting 
polarised Zeeman signatures in spectral lines, are sensi- 
tive to the vector properties of the magnetic field but en- 
counter difficulties in reconstructing highly-complex multi- 
polar field structures. Until now, circular spectropolarimetry 
has been used, for cTTSs, to detect and estimate magnetic 
fields in accretion funnels of cTTSs (|Johns-Krull et al.lll999l : 



IValenti Johns- Krulll 2004 ) . but not the photosphe ric fields 
fe.g. JValenti Johns-KT^l2004l : IPaou et al.ll2006l ). 

Here we report a new study of the cTTS V2129 Oph 
(SR 9, AS 207A, GY 319, ROX 29, HBC 264), the brightest 
TTS in the p Oph cloud, using ESPaDOnS, the new gen- 
eration spectropolarimeter recently installed at the 3.6-m 
Canada-France-Hawaii Telescope (CFHT). After briefly re- 
viewing the fundamental parameters of V2129 Oph (Sec. [2]), 
we present our observations (Sec. [3]) and describe the shape 
and rotational modulation of the Zeeman signatures de- 
tected in both photospheric and emission lines (SecH]). Us- 
ing tomographic techniques, we model the distribution of 
active regions, accretion spots and magnetic fields at the 
surface of V2129 Oph (Sees. [5] & [6]). As an illustration, we 
provide a first attempt at modelling the magnetosphere and 
accretion funnels of V2129 Oph ((SecE)). We finally discuss 
the implications of our results for stellar formation (Sec. [8]). 



2 THE CTTS V2129 OPH 

Several estimates of the photospheric temperature Teg o f 
V21 29 Oph are available in the li terature. Bouvierl (|l990l ) 
and iBouvier Appenzellerl (19921 ) initially suggested that 
Teff - 4200 K. IPadgettI (|l99d ) and [Eisner et al.1 ()2005l ) pro- 
posed hlgh^£_tem2era^ a nd 4400 K respectively ) 
while iGeoffrav MoninI (|200ll ) and'Poppm ann et al.1 (|2003l ) 
concluded that the star is actually cooler (4000 and 3400 K 
respectively). These differences probably arise because the 
spectral energy distribution includes a strong infrared con- 
tribution from the accretion disc (Eisner et al. 2005) and, 
to a lesser extent, from its companion stai0. The photo- 
metric colours of V2129 Oph are, moreover, affected by in- 
terstellar reddening, making photometric temperature esti- 
mates slightly unreliable. The spectroscopic measurements 



by IPadgettI (|l996l ) and iBouvier Appenzellerl (|l992l ) ap- 
pear safest; we thus choose 4500 K as the most likely esti- 
mate, with a typical error bar of about 200 K. 

From this temperature estimate, and given the B — V, 
V — Rc and V — Ic photometric colours of V2129 Op h 
(respectively equal to 1.25, 0.80 and 1.60, iBouvierl [T990l ) . 
we infer from Kurucz synthe tic colours (including interstel- 
lar reddening, iKuruczl 1 19931 ) that the B — V colour ex- 
cess E{B — V) 0.2 (i.e. that the interstellar extinction 
Av — 0.6) and that the bolometric correction (including in- 
terstellar extinction) is —1.3. Although V2129 Oph was too 
faint to be observed with Hipparcos (and thus to have its 
parallax estimated accurately) , its membership of the p Oph 
dark cloud ind icates that its distance is about 140 di 10 pc 
(jBontemps et a l. 2001). 

' Bouvierl (|l990) and ISchevchenko HerbstI (|l998l ) find 
that the maximum brightness of V2129 Oph (at minimum 
spot coverage) is about 11.20; its unspotted V magnitude 
is thus 11.0 =b 0.2, given that active stars usually feature 
as m uch as 20% spot co verage even at brightness maximum 
fe.g. JPonati et aLlliooi ). From this, we obtain that the bolo- 
metric magnitude of V2129 Oph is 4.0 =b 0.25, i.e. that its 
logarithmic luminosity (with respect to the Sun) is 0.3 ±0.1. 
Given the assumed Tgff , we derive that V2129 Oph has a ra- 
dius of = 2.4zb0.3 R©. In Sec. [5] we measure the projected 
equatorial rotation velocity to be vsini = 14.5 ±0.3 kms~^, 
where i is the angle between the rotation axis and the 
Une of sight. This agrees well with published estimate s of 
15.2 ±0.9 kms~^ a nd 15.8±1.5 kms"^ by Eisner et al. 20051 
and[B ouvie 3 ll990l respectively. From v sin i and the rotation 
period of 6.53 d (jSchevchenko k HerbstI 1 19981 ). we obtain 
that R^sini = 1.87 ± 0.04 R©. We infer i 50°, in good 
agreement with the value of 45° derived from tomographic 
modelling of our data (see Sec. [5]). 

Fitting the evolutionary models of ISiess et aD (|2000l ) 
to these parameters, we infer that V2129 Oph is a 1.35 ± 
0.15 Mq star with an age of about 2 Myr and is no 
longer fully convective, with a small radiative core of mass 
^ 0.1 and radius 0.2 R^. 

From the equivalent width of the 866 nm Ca ii line emis- 
sion core (about 0.05 nm, see Sec. 2] implying a line emis- 
sion flux of 10~^ ergs~^cm~^), we obtain that the mass- 
accretion rate M of y2129 Oph is about 4 x 10"^ M© yr"^ 
(iMohantv et al.ll2005l ). lEisner et all (|2005l ) suggests that M 
is significantly larger, about 3 x 10~^ Moyr"^. We assume 
here that M ^ 10~^ M0yr~^, i.e., a conservative compro- 
mise between both estimates and a typical value for a mildly 
accreting cTTS. 

Rotational cycles E are computed according to the 
ephemeris: 



HJD = 2453540.0 + 6.53^, 



(1) 



where the rotation perio d is that determined by 
ISchevchenko HerbstI (|l998l ). 



^ V2129 Oph is a distant binary whose secondary star, pre- 
sumably a very low-ma ss star or a brown dwarf according to 
iGeoffrav MoninI (|200ll) , is about 50 times fainter than the cTTS 
primary in the V band. 



3 OBSERVATIONS 

Spectropolarimetric observations of V2129 Oph were col- 
lected with ESPaDOnS in 2005 June. The ESPaDOnS spec- 
tra span the whole optical domain (from 370 to 1,000 nm) 
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Table 1. Journal of observations. Columns 1-4 sequentially list 
the UT date, the heliocentric Julian date and UT time (both at 
mid-exposure), and the peak signal to noise ratio (per 2.6 kms"-*^ 
velocity bin) of each observation. Column 5 lists the rms noise 
level (relative to the unpolarized continuum level Ic and per 
1.8 kms"-*^ velocity bin) in the circular polarization profile pro- 
duced by Least-Squares Deconvolution (see SecO, while columns 
6 indicates the rotational cycle associated with each exposure (us- 
ing the ephemeris given by Eq[T}. Exposure times are equal to 
4 X 600 s for all observations. 



Date HJD UT S/N ctlsd Cycle 

(2005) (2,453,000+) (h:m:s) (10-^/c) 

Jun 19 540.86507 08:39:04 100 5.5 0.132 



Jun 20 541.80878 07:18:04 100 4.9 0.276 

Jun 20 541.83908 08:01:42 100 4.9 0.281 

Jun 20 541.86941 08:45:23 90 5.8 0.286 



Jun 21 542.81897 07:32:49 110 4.1 0.431 

Jun 21 542.84993 08:17:24 110 4.2 0.436 

Jun 21 542.88152 09:02:54 110 4.1 0.441 



Jun 22 543.89201 09:18:04 70 7.2 0.595 

Jun 22 543.92304 10:02:46 70 6.5 0.600 

Jun 22 543.95467 10:48:19 90 5.3 0.605 



Jun 23 544.87704 08:56:36 120 3.9 0.746 

Jun 23 544.90808 09:41:17 120 3.9 0.751 

Jun 23 544.93912 10:25:60 110 4.1 0.756 



Jun 24 545.88764 09:11:56 120 3.8 0.901 

Jun 24 545.91867 09:56:37 120 3.9 0.906 

Jun 24 545.94970 10:41:19 110 4.2 0.911 



Jun 25 546.87674 08:56:19 120 3.8 1.052 

Jun 25 546.90776 09:40:59 120 3.7 1.057 

Jun 25 546.93878 10:25:40 120 3.8 1.062 



Jun 26 547.88109 09:02:40 120 3.9 1.206 

Jun 26 547.91212 09:47:21 120 3.9 1.211 

Jun 26 547.94317 10:32:03 110 4.0 1.216 



at a resolving power of 65,000. A total of 22 circular- 
polarisation sequences were collected on 8 consecutive 
nights, each sequence consisting of 4 individual subexpo- 
sures taken in different polarimeter configurations. The full 
journal of observations is present ed in Table [T ] The ex- 
traction procedure is described in IPonati et al.l (| 19971 , see 
also Donati et al., 2007, in prep., for further details) and 
is carried out with Libre ESpRIT, a fully automatic re- 
duction package/pipeline installed at CFHT for optimal ex- 
traction of ESPaDOnS unpolarised (Stokes /) and circularly 
polarised (Stokes V) spectra. The peak signal-to-noise ratios 
per 2.6 kms~^ velocity bin are about 1 00 on avera^qj. 

Least-Squares Deconvolution fLSD: lDonati et al.lll997l l 
was applied to all observations. The line list we employed for 
LSD is compu ted from an AtlasQ LTE model atmosphere 
(|Kuruc j \l9m and corresponds to a K5IV spectral type 
(Teff = 4, 500 K and log g = 3.5) appropriate for V2129 Oph 
(see Sec. [2)). We selected only moderate to strong spectral 

At the time of our observations, ESPaDOnS suffered a 1.3-mag 
light loss due to severe damage to the external jacket of optical 
fibres linking the polarimeter with the spectrograph. This has 
since been repaired. 
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LSD profiles, V2129 Oph, ESPaDOnS, 2005 Jun. 24 




Velocity (knn/s) 

Figure 1. LSD circularly-polarized and unpolarized profiles of 
V2129 Oph (top, bottom curves respectively) on 2005 June 24 
(phase 0.905). The mean polarization profile is expanded by a 
factor of 25 and shifted upwards by 1.05 for display purposes. 



Unes whose synthetic profiles had line-to-continuum core de- 
pressions larger than 40% neglecting all non-thermal broad- 
ening mechanisms. We omitted the spectral regions within 
strong Unes not formed mostly in the photosphere, such as 
the Balmer and He lines, and the Ca ii H, K and infrared 
triplet (IRT) lines. Altogether, about 8,500 spectral features 
are used in this process, most of them from Fe L Expressed 
in units of the unpolarised continuum level /c, the average 
noise levels of the resulting LSD signatures range from 3.7 
to 7.2x10"^ per 1.8 kms~^ velocity bin. They are listed in 
Tabled 

Zeeman signatures, featuring full amplitudes of about 
0.5%, are clearly detected in the LSD profiles of all spectra 
(e.g., see Fig. [1]). Circular polarisation is also detected in 
most emission lines, and in particular in the 587.562 nm 
He I line and in the Ca ii (IRT ) known as ^ood trac- 
ers of magnetospheric accr etion ([Johns- Krull et al ] Il999l : 
IValenti Johns-Krulll [20041 ). To increase S/N, LSD pro- 
files from data collected within each night were averaged 
together. Given the very small phase range (typically 1%) 
over which nightly observations were carried out, no blur- 
ring of the Zeeman signatures is expected to result from this 
procedure. The resulting set of LSD Stokes / and V profiles 
(phased according to the ephemeris of eqn. [T]) is shown in 
Fig. [2] while that corresponding to the Ca ii IK^ and He i 
emission lines are shown in Figs.[3]and|4l The corresponding 
longitudinal fields, co mputed from the firs t order moment of 
the Stokes V profiles (|Donati et al.lll99"7h , are listed in Ta- 
ble H 



4 ROTATIONAL MODULATION 

In this section, we examine the temporal variations and 
shapes of the Stokes / and V LSD profiles and emission 



^ Note that the 3 components of the Ca ll IRT were averaged 
together in a single profile to increase S/N further. 
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Figure 2. Average Stokes / (left panel) and Stokes V LSD pro- 
files (right panel) of photospheric lines, for each observing night 
(top to bottom). The mean Stokes / LSD profile (averaged over 
the complete set, thin line) is also shown to emphasise temporal 
variations. The rotational cycle associated with each observation 
is noted next to each profile; 3cr error bars are also shown to the 
left of each Stokes V profile. All profiles are plotted in the velocity 
rest frame of V2129 Oph. 

lines throughout the period of our observations. In partic- 
ular, we show that most of the temporal variations of line 
profiles are attributable to rotational modulation. 



4.1 Photospheric lines and accretion proxies 

The Ca ii IRT and He i emission lines exhibit the strongest 
and simplest evolution over the period of our observations. 
The fluxes and the amplitudes of the Zeeman signatures in 
the emission lines are strongest on nights 5 and 6, and weak- 
est at the beginning and end of our monitoring (see Figs. [3] 
and m right panels) . The typical timescale of this fluctu- 
ation is compatible with the ro t ation period of 6.53 d de- 
rived by ISchevchenko HerbstI (|l998[ ) . Longitudinal fields 
derived from both emission lines (see Table [2]) are plotted 
in Fig. [5] as a function of rotation cycle. These field esti- 
mates return to their initial values (within error bars) once 
the star has completed one full rotational cycle, indicating 
that the temporal fluctuations of the Zeeman signatures are 
mostly attributable to rotational modulation. It also sug- 
gests that the magnetic fields are mostly monopolar over 
the formation region of the Zeeman signatures, which we 
identify with the accretion spots at the footpoints of accre- 
tion funnels. Very similar behaviour is reported for other 
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Figure 3. Same as Fig. [2]for the Ca ii IRT lines. 
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Figure 4. Same as Fig. [2] for the He I (right) lines. 
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Table 2. Longitudinal magnetic field of V2129 Oph, as estimated 
from LSD profiles (^lsd, col. 2), the Call IRT emission core 
(^IRT, col. 3) and the He I line (^He I' Column 7 lists the 

veiling parameter r at each phase, defined as the relative decrease 
in the depth of photospheric lines (with respect to their average 
intensity) and estimated from LSD Stokes / profiles. Rotational 
cycles (col. 1) are computed according to Eq[T] 



Cycle 


^LSD 


^IRT 


^Hei 


r 




(G) 


(kC) 


(kC) 




0.132 


-97± 17 


0.23 ±0.07 


0.47 ±0.62 


±0.01 


0.281 


-6 ±9 


0.24 ±0.05 


0.71 ±0.39 


±0.02 


0.436 


-14 ±8 


0.33 ±0.03 


0.78 ±0.17 


0.00 


0.601 


-35± 11 


0.44 ±0.04 


1.10 ±0.19 


-0.01 


0.751 


+51 ±7 


0.64 ±0.03 


1.11 ±0.09 


±0.02 


0.906 


-94 ±7 


0.67 ±0.02 


1.22 ±0.08 


-0.06 


1.057 


-180 ±7 


0.30 ±0.03 


0.91 ±0.15 


0.00 


1.211 


-39 ±7 


0.28 ±0.03 


0.72 ±0.22 


±0.02 




Rotational cycle 

Figure 5. Longitudinal field estimates derived from the Zeeman 
signatures of the Ca ll IRT (open squares) and He I (filled dots) 
emission lines, as a function of rotational cycle; ±lcr error bars 
are also shown. 



cTTSs fe.g.. lValenti Johns-Kml]|l2004l l. further strength- 
ening our interpretation. 

The equivalent width of the Ca ii IRT line emission 
varies from 14 to 21 kms~^ (0.04 to 0.06 nm), while the He i 
line emission varies from 6 to 18 kms~^ (0.012 to 0.036 nm). 
We note that the variation of this emission strength is ap- 
parently not exactly periodic. While the Ca ii IRT and He i 
emission strengths are equal to 15 and 6 kms~^ respectively 
at the beginning of our observations, they are slightly larger 
(by about 2 and 1 kms~^ respectively) one complete rota- 
tional cycle later. This intrinsic variability of the emission 
strength is however significantly smaller than the rotational 
modulation itself. 

We suggest that the Ca ii IRT and He i emission lines 
comprise two physically distinct components. We attribute 
the first of these, the accretion component, to localised ac- 
cretion spots at the surface of the star whose visibility varies 
as the star rotates, giving rise to rotational modulation of 
the emission. A second, chromospheric component arises in 
basal chromospheric emission distributed more or less evenly 
over the surface of the star, producing a time-independent 
emission component. Assuming that the minimum line emis- 
sion strength over the rotation cycle roughly corresponds to 
the chromospheric component, we obtain that the ratio of 
the accretion component to the basal chromospheric emis- 
sion is about 1:2 and 2:1 for the Ca ii IRT and He i lines 
respectively. If we further assume that the chromospheric 
component is mostly unpolarised while Stokes V signatures 
arise mainly in the accretion component, we deduce that the 
Zeeman signatures are diluted within the full line emission 
fluxes by factors of 1:3 and 2:3 for the Ca ii IRT and He i 
lines respectively. The longitudinal field strengths derived 
from each line need to be scaled up by factors of 3:1 and 3:2 
to recover the longitudinal field corresponding to the accre- 
tion component only. In this way we estimate the projected 
field strength within accretion spots at maximum visibility 
to be about 2 kG. 

The Stokes / and V profiles of the LSD photospheric- 
line profiles also show clear temporal variations (see Fig.[2|). 
These are more complex than those traced by the Ca ii IRT 
and He i emission lines. While the Zeeman signatures of 
the emission lines are mostly monopolar and keep the same 



(positive) polarity throughout the whole rotation cycle, the 
photospheric Stokes V profiles exhibit several sign switches 
across the rotation profile (e.g., up to 3 at phases 0.751 and 
0.905) and change polarity with time (see Table [2])- The 
longitudinal field strengths associated with the photospheric 
Zeeman signatures are usually less than 100 G (except on ro- 
tation cycle 1.057), much weaker than those of the Ca ii IRT 
and He i emission lines. This indicates clearly that the un- 
derlying magnetic topology associated with the photospheric 
Zeeman signatures is tangled and mutipolar, suggesting that 
emission and photospheric lines do not form over the same 
regions of the stellar surface. Given the complexity of the 
photospheric Zeeman signatures, a quick glance at the data 
cannot demonstrate conclusively that the observed modula- 
tion is consistent with pure rotational modulation. We can, 
however, verify that the shape of the Zeeman signature at 
the begining of our observations (phase 0.132 and 0.280) 
is similar to that observed one rotation later (phase 1.211). 
Definite confirmation that this variability is entirely compat- 
ible with rotational modulation needs sophisticated stellar 
surface imaging tools; this is achieved in Sec. [H] 

The photospheric Stokes / LSD profile displays only 
moderate variability, which is confined mainly to the line 
core. The equivalent width of the photospheric line pro- 
file varies by no more than a few percent (see Table [2)) , in 
agreement with previous publishe d reports on V2129 Oph 
(|Padgettl [l996l : lEisner et al"1l2005l ). This suggests that ac- 
cretion spots on V2129 Oph are not hot enough to produce 
significant continuous "veiling" emission at the photospheric 
level (and hence in LSD profiles), a reasonable conclusion for 
a mildly accreting cTTS. The photospheric Stokes / profiles 
exhibit shape variations reminiscent of the distortions in- 
duce d by cool surface spots in spectral lines of active stars 
(e.g., Donati e t al. 2003). Given the small amplitude of these 
distortions and the relatively sparse rotational sampling of 
our data set, we again need stellar surface imaging tools to 
confirm that the profile distortions we detect are compati- 
ble with rotational modulation and attributable to surface 
spots; this is presented in Sec. [S] 

The average radial velocity of the photospheric lines is 
—7.0 kms~^, which we take in the following to be the he- 
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Figure 6. Same as Figs.[2]for Ha (left) and H/3 (right). 
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liocentric radial velocity of the stellar rest frame. The Ca ii 
IRT emission core is centred at —6 kms~^ on average. This 
is only very slightly redshifted with respect to the stellar 
rest frame, and varies by less than ±1 kms~^ over the ro- 
tational cycle (see Fig. [3]). The width of the emission core 
(full width at half maximum of 22 kms~^) is comparable to 
the rotational broadening of the star. 

The He l emission line centroid is more significantly red- 
shifted relative to the stellar rest frame than the Ca ii IRT 
emission core, lying at kms~^ on average. Its full- width at 
half- maximum of 33 kms~^ is also significantly larger than 
that of the Ca ii emission core. Both discrepancies are partly 
attributable to the He i line being composed of 6 different 
transitions, the last of which redshifted from the main group 
by about 20 kms~^. The He i centroid velocity varies in ve- 
locity by up to 2.5 kms~^ with rotation phase, moving from 
blue to red between phase 0.45 and 1.05 and crossing the line 
centre at phase 0.75 (see Fig.|4]). This is greater than the ve- 
locity shifts of the Ca ii IRT emission core. It is qualitatively 
consistent with the simple two-component model of emission 
lines proposed above and suggests that rotational modula- 
tion is mostly caused by the accretion component, which is 
stronger in the He i line than in the Ca ii IRT emission line. 

The widths and redshifts of the Ca ii and He i Zeeman 
signatures are similar to those of the unpolarised emission 
profiles. The Stokes V profiles of the Ca ii IRT emission 
cores are mainly antisymmetric with respect to the line cen- 
tre. The Zeeman signature of the He i line departs signifi- 
cantly from antisymmetry, with a blue lobe both stronger 
and narrower than the red lobe. This suggests that the He i 
line forms in a region featuring non-zero velocity gradients. 



4.2 Balmer lines 

Ha and H/3 lines exhibit strong emission with average equiv- 
alent widths of 900 and 200 kms~\ (i.e., 2 and 0.32 nm) 
respectively. They exhibit clear rotational modulation and 
convey information that Ca ii IRT and He i emission lines 
do not contain (see Fig. [6|). One obvious difference is that 
their unpolarised emission profiles are significantly broader 
(with full widths at half maximum of 190 and 90 kms~^ 
for Ha and H/3 respectively) than their Zeeman signatures 
(70 kms~^ wide). They also include a high- velocity compo- 
nent in their red wings (between +100 and +300 kms""*^) 
that alternately appears in absorption and emission (at cy- 
cles 0.75 and 1.21 respectively) but does not show up in the 
Stokes V profile. Both points are readily visible on the aver- 
age emission and standard deviation profiles of Balmer lines 
(see top panel of Fig. [7] for Ha to H7). 

Autocorrelation matrices of Balmer emission profiles 
(see bottom panel of Fig. [T]) reveal that the high-velocity 
redshifted component is anticorrelated with the central emis- 
sion peak. We propose that the central Balmer emission 
peak, which varies in phase with the Ca 11 and He i emission, 
is partly produced in the accretion shock. From its signifi- 
cant blueshift (7 kms~^ with respect to the stellar velocity 
rest frame for Ha and H/3), the Balmer emission likely in- 
cludes a wind component as well, strongest in Ha. The anti- 
correlated modulation of the high-velocity redshifted com- 
ponent points to a different origin; we suggest it is due to 
free- falling material in the preshock region of the accretion 
funnels. 
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Figure 7. Mean Ho; (left), ¥lf3 (middle) and H7 (right) profiles (top panel, full line) and associated standard deviations (top panel, dashed 
line), with corresponding autocorrelation matrices (bottom panel). Note that the standard deviation profiles are multiplied by 3 and 
shifted upwards by 1 for display purposes. In the lower panels, white indicates perfect correlation and black means perfect anticorrelation. 



5 ACTIVE REGIONS AND ACCRETION 
SPOTS 

From the overview of the collected emission- line data in 
Sec. m we infer that rotational modulation is responsible 
for most of the observed variability in the Ca 11 IRT and 
He I emission lines. We demonstrate in the present and fol- 
lowing sections that this is also the case for the Stokes / 
and V LSD profiles of photospheric lines. As a first step, 
we use Doppler tomography to reconstruct the locations of 
dark spots on the surface of V2129 Oph from the time se- 
ries of unpolarised photospheric line profiles. We also apply 
this Doppler imaging technique to the emission-line profiles, 
using them as accretion proxies to determine where the ac- 
cretion spots are located. 



5.1 Dark spots 

As we noted in Sec. IH the distortions seen in the Stokes 
/ LSD profiles of V2129 Oph resemble those caused by 
cool photospheric feat u res on the surfaces of active stars. 
ISchevchenko HerbstI (|l998l ) reach a similar conclusion 
from multicolour photometry. In the following, we model the 
observed LSD profile variations using tomographic imaging, 
assuming that dark surface features are present at the sur- 
face of V2129 Oph. We thus demonstrate that the temporal 
variations in the LSD profiles are mainly attributable to ro- 
tational modulation, and determine where dark spots are 
located on the surface of V2129 Oph. 

We remove the veiling by simply scaling all profiles 
to the same equivalent width, an d use the maximum en- 
tropy image reconstruction code of lBrown et al and 



iDonati Brownl (|l997l ) to adjust the surface image itera- 
tively in order to optimise the fit to the series of Stokes / 
profiles. The image quantity we reconstruct is the local sur- 
face brightness (relative to the quiet photosphere), varying 
from 1 (no spot) to (no light). The local profile, assumed 
to be constant over the whole star, is obtained by match - 
ing a Unno-Rachkovski profile (|Landi Degl'Innocenti|[20o3 l 
with no magnetic fields to the LSD profile of a K5 template 
star (61 Cyg) computed with the same line list as that used 
for V2129 Oph. The synthetic profiles we finally obtain (see 
Fig. [SI left panel) fit the observed Stokes / LSD profiles at a 
level of about 0.15% or S/N=700 and confirms that the pro- 
file distortions are attributable to cool spots at the surface 
of the star. 

Given the relative sparseness of our data set, one may 
argue that the success at fitting the Stokes / LSD profiles 
is no more than a coincidence and does not truly demon- 
strate that the profiles are rotationally modulated. To check 
this, we varied the assumed rotation period of the star over 
a small interval about the value of Schevchenko & Herbst] 
(| 19981 1 and reconstructed brightness images for each of 
the assumed rotation periods. We anticipate that the re- 
constructed image will show minimal spot coverage at 
the correct rotation period when the data are fitted at 
a fixed x^, (or equivalently will yield the best fit to the 
data for a given spot coverage), if the profile variability 
is truly due to rotational modulation. This technique has 
also been used to derive constraints on the amount of dif- 
fer entia]_^otaAi ion^heaxm 

( Do nati et al.1 |200Q| : IPetit et al.( l2QQ2l : iDonati et all l2003l ). 
Assuming solid body rotation, the optimal period we find 
is 6.56 =b 0.02 d, compatible with the value of 6.53 d derived 
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Figure 8. Observed LSD Stokes / (left panel) and Ca ii IRT 
emission core (right panel) of V2129 Oph (thin line) along with 
the maximum entropy fit to the data (thick line). The rotational 
phase and cycle of each observation is written next to each profile. 



bv lSchevchenko HerbstI (|l998l ). Assume now that the sur- 
face of the star rotates differentially, with the rotation rate 
varying with latitude as Qeq — dQsin^ 0, Qeq being the an- 
gular rotation rate at the equator and dQ the difference in 
angular rotation rate between the equator and pole. We then 
find that optimal fits to the Stokes / data are achieved for 
Qeq = 0.964zb0.003 radd"^ and dQ = 0.054zb0.019 radd"\ 
corresponding to rotation periods of 6.52 and 6.90 d for the 
equator and pole respectively. The error bar on dQ is too 
large to detect unambiguously differential rotation at the 
surface of the star; we can however exclude values of dQ 
that are either smaller than or larger than 0.1 radd~^ 
(i.e., about twice the shear at the surface of the Sun). As 
a by-product, we very clearly demonstrate that the tempo- 
ral variability of Stokes / LSD profiles is due to rotational 
modulation. 

The resulting map (see Fig. (9] left panel) includes one 
cool spot close to the visible pole and several additional ones 
at low to intermediate latitudes. The main polar spot is un- 
doubtedly real. The low- latitude appendages at phases 0.23 
and 0.95 are also well constrained by observed profiles fea- 
turing obvious line asymmetries. The two features centred 
on phase 0.44 and 0.62 are only associated with low level 
line asymmetries and their reality is less well established. 
The total area covered by photospheric cool spots in the re- 
constructed image is about 7%, with about 5% for the main 
polar spot only. We emphasize that this total spot coverage 
should be regarded as a lower limit. Doppler imaging recov- 



ers successfully the most significant brightness features, such 
as the cool polar spot seen here. Even if the stellar surface is 
pepp ered with spo ts much smaller than our resolution limit 
(e.g., Ijeffers et al ]|2006), such features will be suppressed, 
especially in stars with moderate rotation velocities. 

We find that vsini = 14.5 ib 0.3 kms~^ and Vrad — 
—7.0 ±0.2 kms~^ at the time of our observations, in perfect 
agreement with the rec ent published estim ates (15.2 zb 0.9 
and -7.17±0.25 km s~\ [Eisner et aLlliooBI ). Given the spec- 
tral resolution of ESPaDOnS (of about 5 kms~^), the longi- 
tudinal resolution of the imaging process at the stellar equa- 
tor is about 20° or 0.05 rotation cycle. This is consistent with 
the longitudinal extent of the brightness features we recov- 
ered. We find that the information content of the brightness 
image is minimised for an axial inclination angle i — 45°, in 
agreement with expectations from spectrophotometric esti- 
mates (see Sec. [2]). 

5.2 Accretion spots 

We proceed in a similar way to model the fluctuations of the 
Ca II IRT emission profiles and to retrieve information about 
the location of accretion spots at the surface of V2129 Oph. 
Thanks to their much smaller noise level, narrower width 
and simpler Zeeman signatures, and despite their weaker 
amount of rotational modulation, the Ca il emission lines 
are a better choice than the He i line for this modelling task. 
The modelling accuracy is, however, limited by the intrinsic 
variability of the unpolarised emission lines. For this mod- 
elling, we use the simple two-component model described 
in Sec. 14] involving a basal chromospheric emission compo- 
nent evenly distributed over the surface of V2129 Oph and 
an accretion emission component arising in local accretion 
spots. 

We now reconstruct an image of the local Ca ii excess 
emission from accretion spots relative to the basal chromo- 
spheric emission. We assume that the Ca ii minimum emis- 
sion profile (at phase 0.28, see Fig. [3]) provides a rough esti- 
mate of the basal chromospheric component. To model this 
chromospheric component, we use a Gaussian local profile 
with a full width at half maximum of 12 kms~^ and an 
equivalent width of 15 kms~^; given the vsini estimate de- 
rived from Stokes / LSD profiles, this model provides a rea- 
sonable match to the minimum emission profile at phase 
0.28. We use the same local profile to model the accretion 
component. The fit we obtain is shown in Fig. [8] (right panel) 
while the reconstructed map is presented in Fig. [9] (right 
panel). As a result of intrinsic variability, the quality of the 
fit is poorer than that for photospheric lines. Using the last 
6 profiles only improves the fit slightly but induces minimal 
changes in the resulting map. We conclude that the recov- 
ered image is robust despite the low-level intrinsic variabil- 
ity. 

The map of Ca ii excess emission at the surface of 
V2129 Oph shares similarities with the distribution of dark 
spots. It shows in particular one main feature close to the 
pole and extending towards lower latitudes between phase 
0.80 and 0.95, directly reflecting (i) the increased Ca ii emis- 
sion at these phases and (ii) the fairly weak velocity change 
of the line emission core. Excess emission from this polar 
accretion spot reaches about 1.5 times that of the average 
quiet chromosphere; its relative area is about 5% of the to- 
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Figure 9. Surface maps of the local surface brightness relative to that of the quiet photosphere (left panel) and the local Ca ll excess 
emission relative to the basal chromospheric emission (right panel) on V2129 Oph, reconstructed from our series of ESPaDOnS LSD 
Stokes / and Can IRT emission line profiles (Fig. [8]). The star is shown in fiattened polar projection down to latitudes of —30°, with 
the equator depicted as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate the phases of observation. 



tal stellar surface. Given the finite resolution of the imaging 
process, the estimate we derive for the fractional area cov- 
ered by accretion spots is likely an upper limit only. 



6 THE LARGE-SCALE MAGNETIC 
TOPOLOGY 

In this section, we use stellar surface imaging to model the 
Zeeman signatures from both LSD profiles and accretion 
proxies. As in Sec. O our aim is to demonstrate that the 
variability of Stokes V signatures is due to rotational mod- 
ulation and to obtain a consistent description of the large- 
scale field topology of V2129 Oph. 

6.1 Model description 

The obvious difference between the sets of Zeeman signa- 
tures from the photospheric LSD profiles and the accretion 
proxies (see Sec. |4]) led us to suggest that photospheric and 
emission lines do not form over the same regions of the stellar 
surface. This is compatible with our findings of Sec. [5] that 
the accretion emission is concentrated in the dark polar re- 
gions. The low surface brightness is expected to suppress the 
photospheric Zeeman signature, rendering the photospheric 
lines insensitive to the field. In the following, we assume this 
to be the case and demonstrate that both sets of Zeeman sig- 
natures are nonetheless consistent with a single large-scale 
magnetic topology. 

As in SecO we use the Ca ii IRT lines as our preferred 
accretion proxy. We describe it with the two-component 
model of Sec. |4l combining a basal chromospheric emission 
component evenly distributed over the star with an addi- 
tional emission component concentrated in local accretion 
spots. The emission-line Zeeman signatures are assumed to 
be associated with the accretion component only, as sup- 
ported by the conclusions of Sec. ID 

The model we adopt involves a vector magnetic field B 



and a local accretion filling-factor / which describes the lo- 
cal sensitivity to accretion proxies and photospheric linefl 
For / = 0, the local area on the protostar's surface con- 
tributes fully to photospheric lines and generates no Ca ii 
excess emission (and only unpolarised chromospheric Ca ii 
emission); for / = 1, the local area does not contribute at 
all to photospheric lines and produces the maximum amount 
of excess Ca ii emission and polarisation (in addition to the 
unpolarised chromospheric Ca ii emission) . Spectral contri- 
butions for intermediate values of / are derived through lin- 
ear combinations between the / = and / = 1 cases. We 
obtain both B and / by fitting the corresponding synthetic 
Stokes V profiles to the observed Zeeman signatures from 
both LSD profiles and Ca ii emission cores. We can also in- 
corporate the additional constraint of fitting the unpolarised 
profiles of the Ca ii emission cores. 

The code we use for fitting B and / is adapted from 
the st ellar surface magn e tic im aging code of iDonati et al.l 
(|200lh and iDonati et al.1 (|2006l ). In this model, the surface 
of the star is divided into a grid of thousands of small sur- 
face pixels, on which we reconstruct B as a 3 component 
vector field and / as a simple scalar field. The magnetic vec- 
tor field is described and computed as a spherical-harmonic 
expansion, whose coefficients are determined through fit- 
ting the simulated Stokes V profiles to the observed data 
sets. More specifically, the field is expressed as the sum of a 
poloidal field and a toroidal field, with the three field com- 
ponents in spherical coordinates being defined through three 
sets of complex spherical harmonics coefficients, a£,m, Pi,m 
and 7£,m- Here i and m denote the order and degree of the 
spherical-harmonic mod^fl- 

For a given magnetic topology, the unpolarised and po- 



Note that the local accretion filling-factor we define here is 
different from the usual 'accretion filling factor' of the cTTS lit- 
erature, i.e., the relative area of the total stellar surface covered 
by accretion spots. 

^ While ai^rn describes the radial field, f^i^rn describes the non- 
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Figure 10. Stokes V Zeeman signatures (thin line) from LSD pro- 
files of photospheric lines (left column) and Ca ll emission cores 
(right column) along with the maximum entropy fit to the data 
(thick line). The rotational phase and cycle of each observation is 
written next to each profile. A 3a error bar is also shown to the 
left of each profile. 



larised spectral contributions to photospheric lines are ob- 
tained through the Unno-Rachkovsky model used in Sec. O 
and assuming a Lande factor and mean wavelength of 1.2 
and 620 nm respectively. For the excess Ca ii emission pro- 
duced in the accretion region, we assume the same Gaus- 
sian model as that of Sec. O with an equivalent width set to 
twice that of the quiet chromosphere profile (thus defining 
the highest possible amount of local excess emission), a unit 
Lande factor and a mean wavelength of 850 nm. 

Given the spatial resolution of the imaging process (see 
Sec. [5]), we find it sufficient to limit the spherical-harmonic 
expansion of the magnetic field components to terms with 
^ < 20 (230 modes altogether). In practice, little improve- 
ment is obtained when adding spherical harmonic terms cor- 
responding to ^ > 15. 



6.2 Modelling results 

We find that the Stokes V signatures of both the photo- 
spheric LSD profiles and the Ca II emission cores can be fit- 
ted simultaneously down to the noise level using the model 
described in Sec. 110! above. Repeating the experiment for a 



radial poloidal field and 7£,rn the toroidal field. More details about 
this description can be found in lDonati et al ] (l2006l V 
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Figure 11. Map of the local accretion filling factors at the surface 
of V2129 Oph, as derived from simultaneously fitting the two 
series of Zeeman signatures (LSD profiles of photospheric lines 
and Ca ll emission cores, see Sec.|4} and assuming that each series 
senses a different spatial region of the stellar surface (see Sec. [6]). 



range of rotation periods (as we did in Sec. [5] for Stokes 
/ profiles only) and assuming solid body rotation yields 
optimal fits to the Stokes V data at a period of 6.62 di 
0.07 d, again consisten t with the 6.53 d period derived by 
ISchevchenko HerbstI (| 19981 ). Assuming now that the star 
rotates differentially, we find that optimal fits to the Stokes 
V profiles are achieved for Qeq = 0.928 it 0.020 radd~^ and 
dO, — —0.038 =b 0.030 radd~^, compatible within 3cr with 
the values derived from the Stokes / profiles. This confirms 
our earlier conclusion from inspection of the profiles in Sec.|4] 
that rotational modulation dominates the temporal variabil- 
ity of the Zeeman signatures. 

The maps of local accretion filling-factors and magnetic 
field derived from the Stokes V data are shown in Figs. [11] 
and [12] (left panel) . They confirm and strengthen the con- 
clusion that accretion is very localised at the stellar surface, 
being concentrated in a high-latitude region where the ra- 
dial field reaches a peak value of 2.1 kC, in agreement with 
the estimate derived in Sec. [4] This accretion region occu- 
pies about 5 percent of the total stellar surface. Including 
a simultaneous fit to the unpolarised Ca ii emission profiles 
has little effect on the result, provided we assume the local 
Ca II excess emission in the accretion regions to be about 
twice the emission from the quiet chromosphere. Note that 
this fractional area is only an upper limit. A similar fit to 
the data can be obtained by assuming that the local Ca ii 
excess emission in the accretion spot is larger than that con- 
sidered here. In this case, we obtain a smaller accretion hot 
spot whose average location is unchanged. The spatial res- 
olution is however high enough to claim that the accretion 
spot cannot be larger than the one we reconstruct. 

As a comparison, we also show in Fig. [12] (right panel) 
the magnetic field topology reconstructed from using the 
photospheric Stokes V LSD profiles alone. Unsurprisingly, 
this topology is essentially the same as that in the left-hand 
panel for all regions where no accretion occurs and for which 
all available information comes from the photospheric lines. 
This map, however, misses most of the magnetic ffux from 
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Figure 12. Reconstructed magnetic topology of V2129 Oph, simultaneously derived from our two series of Zeeman signatures (LSD 
profiles of photospheric lines and Ca ll emission cores, see Sec. jj] left panel) and when using LSD Stokes V profiles from photospheric 
lines only (right panel). The image on the right misses most of the 2 kG high-latitude region of positive radial field, in which accretion 
occurs. The three components of the field in spherical coordinates are displayed (from top to bottom) in both cases, with field fiuxes 
labelled in G. 
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Figure 13. Decomposition of the azimuthal and meridional components of the reconstructed field (see Fig. [TJl left column) into their 
non-radial poloidal (left column) and toroidal (right column) terms, corresponding respectively to spherical-harmonic coefficients Pi^rn 
and 7^,^. 



the high- latitude accretion spot, which now shows up as a 
weak positive field feature only. 

6.3 Field topology 

Even though it is dominated by a 2 kG high-latitude posi- 
tive radial field feature, the magnetic topology nevertheless 
differs significantly from a dipole and exhibits a number of 
clear specific features. In particular, we find that the fit to 
the data (and in particular to the Stokes V LSD profiles) 
is significantly poorer if we assume that the field is purely 
poloidal, implying that the toroidal field component we re- 
construct is likely to be real. The spatial structure of this 
toroidal component is shown on Fig. 1131 fright column) along 
with the non-radial poloidal field component (left column). 

The large-scale field topology includes a strong axisym- 
metric component. The poloidal field component features a 
strong positive radial field spot near the pole, surrounded by 
a negative radial field ring located slightly above the equa- 



tor (see top left image of Fig. [12]), suggesting a dominant 
axisymmetric mode with i higher than 1 . This is also visible 
from the dual ring structure in the meridional component 
of the poloidal field (see bottom left image of Fig. [13]). The 
toroidal field topology is simple (see right panel of Fig. [T3)) 
and comprises a ring of counter-clockwise oriented field. This 
unipolar ring is tilted with respect to the rotation axis to- 
wards phase 0.55 and passes through the rotation pole, giv- 
ing it a complex multipolar appearance when plotted in 
spherical coordinates. 

Since the lower hemisphere of V2129 Oph is poorly 
constrained by observations (southern surface features be- 
ing strongly limb-darkened and only visible for a short frac- 
tion of the rotation cycle) , there is no direct way of working 
out whether the large-scale field topology is mainly symmet- 
ric or antisymmetric with respect to the centre of the star, 
or whether it features a mixed combination of symmetric 
and antisymmetric modes. By forcing the code towards ei- 
ther symmetric or antisymmetric topologies or by letting it 
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Figure 14. Modulus (in G) of the spherical-harmonic complex coefficients a^^rn: Pi,m and 7£,rn (from left to right, see text) for the 
reconstructed field topology of V2129 Oph (Fig. [TJl left column), as a function of mode degree £ and order m. In this case, the imaging 
code was pushed towards antisymmetric field configurations by favouring modes with odd i values. Modes with ^ < 10 only are displayed 
here. 



choose by itself, we find that all three options are possible 
and correspond to magnetic topologies that are very simi- 
lar over the visible hemisphere and that fit the observations 
equally well. While the mainly symmetric and antisymmet- 
ric magnetic topologies we reconstruct both include a strong 
field region near the invisible pole (mirroring that detected 
close to the visible pole, with equal and opposite polarities 
respectively), the third solution mixes both symmetric and 
antisymmetric solutions in roughly equal proportions and 
produces a topology with no strong field region in the in- 
visible hemisphere. Assuming that magnetic topologies of 
cTTSs are either mainly symmetric or antisymmetric (most 
of the m featuring strong unipolar field regions close to th e 
pole, V alenti Johns-Krullll20Q4l : ISvmington etHI l2005al ) . 
we find that the mainly antisymmetric field topology is more 
likely as it requires a field with lower contrast for the same 
quality of the fit to the data. 

In this antisymmetric case, the dominant mode of the 
reconstructed large-scale field corresponds to an octupole 
(^ = 3), with only little energy reconstructed in the dipole 
field configuration. Fig. [TJ] shows how the reconstructed 
magnetic energy is distributed among the various modes. 
While the toroidal field term is mainly confined to low order 
modes (mostly ^ = 1), the poloidal field extends to modes 
with at least i — 5 - with a maximum at ^ = 3 - and tends 
to concentrate mostly in axisymmetric modes with small m 
values. Note that, given the limited resolution of our imaging 
process, we miss all magnetic flux at small spatial scales (e.g. 
in the form of closed bipolar groups) that may be present 
on V2129 Oph. 



7 MAGNETOSPHERIC ACCRETION AND 
CORONA 

To illustrate how accretion proceeds between the inner disc 
and the surface of the star, we extrapolate the reconstructed 
magnetic fleld over the whole magnet osphere. To do this, we 
assume that the 3D fleld topology is mainly potential and 
becomes radial beyond a certain magnetospheric radius ta 
from the star, to mimic the opening of the largest mag- 



netic loops under the co ronal pressure (|jardine et al.|[200^ , 
l2006l : iGregorv et aD[2006l ). In non-accreting stars, this dis- 
tance is usually assumed to be smaller than or equal to the 
corotation radius (rc) at which the Keplerian orbital pe- 
riod equals the stellar rotation period. In cTTSs however, 
the magnetic fleld of the protostar is presumably clearing 
out the central part of the accretion disc, suggesting that 
the magnetosphere sho uld extend as far as the inner disc 
rim ('jardi ne et al.l l2006) . The derived magnetospheric maps 
(see Fig. [15] for two possible values of ta) clearly show that 
the fleld topology is complex, with an intricate network of 
closed loops at the surface of the star. As ta is increased, 
the magnetosphere is dominated by more extended open and 
closed fleld lines. 

For several reasons, our description of the magneto- 
sphere is only an approximation. The innermost fleld ge- 
ometry is likely to be more complex than our reconstruc- 
tion, as small-scale magnetic structures at the surface of the 
star remain undetected through spectropolarimetry. This is 
not, however, a crucial issue for the present study, which 
focusses mainly on the medium- and large-scale magnetic 
topology of V2129 Oph. Moreover, the true magnetic topol- 
ogy of V2129 Oph is likely not potential given the strong 
plasma flows linking the disc to the stellar surface and the 
signiflcant toroidal fleld component detected at the surface 
of V2129 Oph. Finally, the magnetic fleld in the disc is also 
expected to interact with t he magnetosphere , at least i n 
its outer regions (e.g.,|von Rekowski Brandenburdl2004l l. 
These complications are not, however, included in most ex- 
isting theoretic al studies on magnetospheric accretion (e.g., 
iRomanova et al. 2003, 2004). We thus believe that, despite 
its limitations, our imaging study should be robust enough 
to yield the rough locations of accretion funnels and to pro- 
vide a useful illustration of the constraints that spectropo- 
larimetric data sets can yield. More quantitative analysis 
and simulations are postponed to forthcoming papers. 

We can estimate where the accretion funnels are located 
and where they are anchored at the surface of the star, by 
identifying those magnetospheric fleld lines that are able to 
accrete material from the disc. Such fleld lines must link the 
star to the disc and intersect the rotational equator with 
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Figure 15. Magnetospheric topology of V2129 Oph as derived from potential extrapolations of the reconstructed surface magnetic 
topology (left panel of Fig. I12|) . The magnetosphere is assumed to extend up to the inner disc radius, equal to 3.4 and 6.8 R~k in the left 
and right panels respectively. The complex magnetic topology close to the surface of the star is very obvious. In both cases, the star is 
shown at rotational phase 0.8. The colour patches at the surface of the star represent the radial component of the field (with red and 
blue corresponding to positive and negative polarities); open and closed field lines are shown in blue and white respectively. 



effective gravity pointing inwards in the co-rotating frame 
of reference (e.g. Gregory et al. 2006, 2007). 

If we assume that ta = 3.4 i?^, we find that disc mate- 
rial accretes mainly on to the radial field spots located close 
to the equator, with only a small fraction of the accretion 
flow reaching the high-latitude spots. This is obviously in 
strong contradiction with the conclusions of Sec. [H] If we 
assume a more extended magnetosphere with ta = 6.8 i?^, 
most of the accreting material flows onto the high-latitude 
accretion spots. Pushing ta further out (to 7.5 i?*) focusses 
all the accreting material at high latitudes only, in much 
better agreement with observations (see Figs [T6l and fTT)) . We 
therefore conclude that ta is closer to 7 i?^ in V2129 Oph 
and likely coincides with rc, equal to 0.075 =b 0.003 AU or 
6.8 ± 1.2 R^. 

In the magnetospheric topology we derive, the northern 
accretion funnel passes in front of the high-latitude accretion 
hot spot between rotation phases 0.7 and 0.8. We suggest 
that free-falling material in this accretion funnel is responsi- 
ble for the high- velocity absorption component observed in 
the red wing of Balmer lines at cycle 0.75 (see Sec. |4|, just 
as solar prominences transiting the solar disc are viewed as 
dark filaments. Similarly, we propose that free-falling plasma 
along the southern accretion funnel that plunges into the in- 
visible accretion spot produces, through scattering, the red- 
shifted emission component that Balmer lines feature at cy- 
cle 1.21. The orientation of the southern accretion funnel is 
indeed ideally oriented for this purpose at this particular ro- 
tation phase. Detailed radiative transfer compu tations (e.g., 
ISvmington et al.l l2005bl : iKurosawa et al.l l2006l ) are needed 
to confirm this. 

The magnetospheric field geometry derived from the re- 



constructed field maps can also be used to model the corona 
of V2129 Oph and estimate the associated X-ray flux. As- 
suming — 1 and an isothermal corona at a temper- 
ature of 20 MK (typical for the hot component of X-ray 
emission fro m^ cTTSs) filled with plasma in hydrostatic eg ui- 
librium (e.g. "jardi ne et al.l 12006: "Gregor v et al.l 12006'). we 
obtain an X-ray emission measure of 4 x 10^^ cm~^ and an 
average (i.e., emission- measure- weighted) coronal density of 
5 X 10^ cm~^. This value is typical for solar-mass cTTSs 
such as V2129 Oph according to the statistics derived from 
the COUP (Chandr a Orion Ultradeep Project) survey (e.g., 
I Jardine et aDl2006l l. It translates into a X-ray luminosity 
of 4 X 10^° ergs~ \ matching well the published value of 
2.5 X 10^° ergs"^ (| Casanova et al.lll995l l. 



8 DISCUSSION 

The spectropolarimetric data we have collected on 
V2129 Oph yields the first realistic model of the large- 
scale magnetic topology on a mildly accreting cTTS. Al- 
though this first study of its kind addresses only one star, 
V2129 Oph can nevertheless be regarded as prototypical of 
mildly accreting cTTSs, exhibiting no behaviour that would 
mark this star out as peculiar. We can thus use it to read- 
dress a number of issues regarding cTTSs. 

8.1 Cool spots and hot accretion regions 

We observe that profile variations due to intrinsic variability 
(e.g., resulting from unsteady accretion flows) are smaller 
than those caused by rotational modulation, for both the 
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Figure 16. Locations of field lines capable of accreting material from the accretion disc to the surface of V2129 Oph, as derived from 
the magnetospheric topologies shown in Fig. [15] Open field lines connecting to the disc are shown in blue, while closed field lines appear 
in white. As in Fig. [15] the magnetosphere is assumed to extend to either 3.4 or 6.8 (left and right panels respectively). The topology 
of the accretion columns becomes much more complex when the magnetosphere is assumed to be smaller. 



photospheric lines and the accretion proxies. This conclusion 
is already robust and confirms previous published results 
based on data sets with similar time samp ling and coverage 
to ours fe.g. JValenti &: Johns-Knil]|2004h . Future observa- 
tions spanning timescales of at least 2 complete rotation 
cycles, would be welcome to establish whether this situation 
is typical. 

The photospheric line profiles of V2129 Oph are vari- 
able in shape and to a lesser extent in strength. The profile 
distortions consist of bumps (rather than dips) and can be 
attributed to the presence of cool spots at the surface of 
the star, rotating in and o ut of the observer's view. Thi s 
confirms the conclusions of ISchevchenko k HerbstI (|l998h . 
who deduced from the photometric colour variations that 
the brightness inhomogeneities on V2129 Oph are mostly 
due to cool and dark spots rather than to bright and hot re- 
gions. With Doppler imaging, we successfully reproduce the 
rotational modulation of photospheric line shapes and find 
that the surface spot distribution is dominated by a cool 
region over the visible pole. In this respect, this spot dis- 
tribution differs from those of cool active stars with similar 
characteristics, which usually do not show polar features at 
rotation periods longer than abou t one week (e.^ ., the Kl 
evolved RS CVn subgiant II Peg, P etit et al.ll2Q06h . 

Accretion spots are also present at the surface of the 
star. Their spectral signatures are not seen at photospheric 
level, but show up clearly in emission proxies such as the 
Ca II IRT and He i lines which usually trace heating of the 
star's upper atmospheric layers. 

These emission proxies include two components: a sta- 
tionary emission component presumed to originate in a 
quiet chromosphere covering the whole stellar surface, and a 



rotationally-modulated emission feature presumably formed 
in local accretion spots. The accretion component exhibits 
only small-amplitude radial velocity modulation, indicating 
that accretion spots are located at high latitudes. This is 
confirmed through tomographic imaging, which shows that 
accretion on the visible hemisphere of V2129 Oph is con- 
centrated in a single high-latitude spot covering less than 
5% of the stellar surface, located within the main cool po- 
lar spot detected at photospheric level. It suggests that the 
heat produced in the accretion shock is not transferred to 
the photosphere efficiently enough to warm it up above the 
temperature of the non-accreting photosphere. We speculate 
that a second accretion spot is present near the invisible 
pole, at a location symmetrical with respect to the centre of 
the star. 

We also find that differential rotation of V2129 Oph is 
comparable with or less tha n that of the Sun, in agreement 
with previous o bservations ([johns-Krulllll996h and theoret- 
ical predictions (jKiiker Rudigei]ll997l ). Our result invali- 
dates in particular early suggestions that cTTSs might be 
strong differential rotators (fSmith..l994 V 



8.2 Magnetic topology and origin of the field 

We observe that rotational modulation dominates the tem- 
poral evolution of the Stokes V profiles of both photospheric 
lines and accretion proxies. Comparing Zeeman signatures 
from both data sets suggests that they form in different 
complementary regions of the stellar surface, with accretion 
proxies tracing magnetic fields at the footpoints of accretion 
funnels and photospheric lines diagnosing fields in the non- 
accreting photosphere. With this assumption, we are able 
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Figure 17. Location of the footpoints of accreting field lines at 
the surface of the star (white circles), assuming that the mag- 
netosphere extends to 6.8 R^,. In this case, accreting field lines 
concentrate almost exclusively in the strong unipolar radial field 
regions at high latitudes. The colours at the surface of the star de- 
pict the radial field component (with red and blue corresponding 
to positive and negative polarities). Note that phase runs back- 
wards with longitude, i.e., decreases from 1 to while longitude 
increases from 0° to 360°. 



to reconstruct, through tomographic imaging, a consistent 
magnetic topology that simultaneously reproduces both sets 
of Zeeman signatures. 

Despite being anchored in 2 kG high-latitude radial 
field spots, the reconstructed magnetic field is more complex 
than the simple dipole assumed in most theoretical mod- 
els of magnetospheric accretion fe.g. jRomanova et al.ll2003l . 
I2OO4I ) . The negative radial field ring encircling the main 
positive radial field pole provides direct evidence that the 
large-scale poloidal field structure corresponds to spherical 
harmonics order ^ > 2. We find that the field of V2129 Oph 
is more likely to be mainly antisymmetric (rather than sym- 
metric) with respect to the centre of the star; the respective 
contributions to the 2 kG radial field region from the ^ = 1 
(tilted dipole), i = 3 (tilted octupole) and ^ = 5 terms are 
about 0.35, 1.2 and 0.6 kG. The reconstructed field also in- 
cludes a toroidal field component in the form of a tilted ring 
of counterclockwise azimuthal field. 

Given that fossil magnetic fields are presumably dis- 
sipated very quickly i n largely convective stars (e.g., 
IChabrier k Kiikerl l2006l \ the most natural assumption is 
that large-scale fields of cTTSs are generated through dy- 
namo action. The mainly axisymmetric field reconstructed 
on V2129 Oph argues in this direction. The magnetic topol- 
ogy seen here, however, differs in important respects from 
those of moderately ro tating cool activ e stars with simi- 
lar characteristics fe.g. JPetit et al.|[2006h . The strong high- 
latitude unipolar field region detected on V2129 Oph and 
the relatively strong i > 3 orders of the large-scale field are 
the most obvious. These differences suggest that cTTSs may 
feature a composite field, including both fossil and dynamo 
components. More data are needed to confirm this and to 
explore the potential impact on stellar formation. 



8.3 Disc-star magnetic coupling 



Several theoretical papers (e^ iKonigll I1991I : [Sh u et al. 
1 19941 : ICameron &; Campbelll Il993l ) studied how the stellar 



magnetic field interacts with the surrounding accretion disc 
and disrupts its vertical structure close to the star. They 
showed further that the balance between accretion torques 
and angular momentum losses causes the rotation of the star 
to evolve towards an equilibrium in which the disc disrup- 
tion radius lies close to ta and just inside the co-rotation 
radius rc- Assuming different models of the star/disc in- 
teraction, these authors derived convenient expressions for 
ta; they proposed that this coupling causes cTTSs to slow 
down to the Keplerian orbital period at ta (thus equating 
rA and rc), explaining why cTTSs are on average rotating 
more slowly than their disc-less equivalents. This scenario is 
refered to as 'disc-locking' in the literature. 

These papers can be used to estimate ta as a function 
of the mass, radius, mass-accretion rate and magnetic dipole 
strength of the protostar. Using values listed in Sec. [2] and 
a dipole strength of 0.35 kG (see Sec. [6l) , we find that ta 
can be as low as 2 R^, if using the model o f iKonigll (|l99lh or 
IShu et al.l (ll994l). i.e., far smaller than rc (about 7 Ri,). The 
model of Cam eron &: Campbelll (|l993h , yielding ta values 
of 5.1 to 8.1 Ri, (in better agreement with rc), appears to 
provide a better quantitative description of the star/disc 
magnetic interaction assuming this coupling is responsible 
for the slow rotation of V2129 Oph. 

The disc-locking scenario was recently criticised on 
various grounds. Observations revealed that surface mag- 
netic strengths on cTTSs (measured from Zeeman mag- 
netic broadening of unpolarised profiles) do not correlate 
well with those expected from theoretical models, cast- 
ing d oubts on the validity of the disc-locki ng mechanism 
(e.g., iBouvier et al.1 l2007l : Ijohns-Krull l2007l ). On theoreti- 
cal grounds, authors claimed that a wind fro m the proto- 
star would blow the field open at 3 i?^ (e.g., ISafie rl 119981 : 
Matt k Pudrit3 l2004l ) and prevent magnetic coupling be- 
tween the star and disc. 

Our study however demonstrates that the most rele- 
vant physical quantity for studying this issue is not the sur- 
face magnetic field itself, but rather the large-scale dipole 
strength on which magnetic coupling mostly depends. On 
V2129 Oph at least, the large- scale dipole component iscom- 
patible with the predictions of" Cameron & CampbelJ (|l993[ ). 
Moreover, our observations establish that accretion occurs 
onto the polar regions of V2129 Oph rather than the equator, 
a phenomenon which we can only explain if the high-latitude 
field links to the accretion disc at a distance of about 7 i?^. 
This is obviously not compatible with the field being blown 
open by the wind at 3 i^^. 



9 CONCLUSION 

In this paper, we report the discovery of magnetic fields 
on the cTTS V2129 Oph using ESPaDOnS, the new stel- 
lar high-resolution spectropolarimeter recently installed at 
CFHT. Circular polarisation Zeeman signatures are de- 
tected in photospheric lines and in emission lines tracing 
magnetospheric accretion. We demonstrate that the tempo- 
ral variations of both unpolarised and circularly polarised 
profiles, monitored over the complete rotational cycle of 
V2129 Oph, are mostly attributable to rotational modula- 
tion. 

From our sets of Zeeman signatures from photospheric 
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lines and accretion proxies simultaneously, we recover the 
medium- and large-scale magnetic topology at the surface 
of V2129 Oph using tomographic imaging tools. In this pro- 
cess, we also successfully derive the location of the accretion 
regions, tracing the footpoints of accretion funnels at the 
surface of the star. We find that the magnetic topology of 
V2129 Oph is significantly more complex than a dipole and 
is dominated by a 1.2 kG octupole, tilted by about 20° with 
respect to the rotation axis. The large-scale dipole compo- 
nent is much smaller, with a polar strengh of only 0.35 kG. 
The magnetic field is mainly antisymmetric with respect to 
the centre of the star, and the accretion regions are found 
to coincide with the two main high-latitude magnetic poles 
(each covering about 5% of the total stellar surface). The 
high-latitude accretion spots apparently coincide with dark 
polar features at photospheric level. The magnetic topology 
of V2129 Oph is unusually complex compared to those of 
non-accreting cool active subgiants with moderate rotation 
periods. 

As an illustration, we also provide a first attempt at 
modelling the magnetospheric topology and accretion fun- 
nels of V2129 Oph using field extrapolation. We find that 
the magnetosphere must extend to distances of about 7 R^, 
to produce accretion funnels that match the observations 
and anchor in the high-latitude accretion spots identified 
at the stellar surface. This distance is roughly equal to the 
corotation radius and n iatche s the theoretical predictions of 
'Cameron & Campbell' ('1993'), suggesting that the star/disc 
magnetic coupUng can possibly explain the slow rotation of 
V2129 Oph. It suggests that magnetic field lines from the 
protostar are apparently capable of coupling to the accretion 
disc beyond 3 Ri, . The accretion geometry we derive is qual- 
itatively consistent with the modulation of Balmer lines and 
consistent with the X-ray coronal fluxes typical of cTTSs. 

More similar data sets are required to check how the 
conclusions of this prototypical study apply to other cTTSs, 
with different masses, rotation and accretion rates. 

ACKNOWLEDGEMENTS 

We thank the CFHT staff for their help during the various 
runs with ESPaDOnS. 

REFERENCES 

Bontemps S., Andre P., Kaas A., Nordh L., Olofsson G., 
Huldtgren M., Abergel A., Blommaert J., Burgdorf M., 
Cesarsky C, Cesarsky C, Copet E., Davies J., Falgarone 
E., Lagache G., Montmerle T., Perault M., 2001, A&A, 
372, 173 

Bouvier J., 1990, A J, 99, 946 

Bouvier J., Alencar S., Harries T., Johns-Krull C, Ro- 
manova M., 2007, in Reipurth B., Jewitt D., Keil K., eds, 
Protostars and Planets V Vol. 951. University of Arizona 
Press, p. 479 
Bouvier J., Appenzeller I., 1992, A&AS, 92, 481 
Brown S., Donati J.-F., Rees D., Semel M., 1991, A&A, 
250, 463 

Cameron A., Campbeh C, 1993, A&A, 274, 309 



Casanova S., Montmerle T., Feigelson E., Andre P., 1995, 

ApJ, 439, 752 
Chabrier G., Kiiker M., 2006, A&A, 446, 1027 
Daou A., Johns-Kruh C, Valenti J., 2006, ApJ, 635, 466 
Donati J.-F., Brown S., 1997, A&A, 326, 1135 
Donati J.-F., Cameron A., Petit P., 2003, MNRAS, 345, 

1187 

Donati J.-F., Cameron A., Semel M., Hussain G., Petit 
P., Carter B., Marsden S., Mengel M., Lopez Ariste A., 
Jeffers S., Rees D., 2003, MNRAS, 345, 1145 

Donati J.-F., Howarth I., Jardine M., Petit P., Catala 

C, Landstreet J., Bouret J.-C, Alecian E., Barnes J., 
Forveille T., Paletou F., Manset N., 2006, MNRAS, 370, 
629 

Donati J.-F., Mengel M., Carter B., Marsden S., Cameron 

A., Wichmann R., 2000, MNRAS, 316, 699 
Donati J.-F., Semel M., Carter B., Rees D., Collier 

Cameron A., 1997, MNRAS, 291, 658 
Donati J.-F., Wade G., Babel J., Henrichs H., de Jong J., 

Harries T., 2001, MNRAS, 326, 1256 
Doppmann G., Jaffe D., White R., 2003, AJ, 126, 3043 
Eisner J., Hillenbrandt L., White R., Akeson R., Sargent 

A., 2005, ApJ, 623, 952 
Geoffray H., Monin J.-L., 2001, A&A, 369, 239 
Gregory S., Jardine M., Simpson I., Donati J.-F., 2006, 

MNRAS, 371, 999 
Gregory S., Jardine M., Wood K., 2007, MNRAS, in press 

(astro-ph/0704.2958) 
Jardine M., Cameron A., Donati J.-F., Gregory S., Wood 

K., 2006, MNRAS, 367, 917 
Jardine M., Wood K., Cameron A., Donati J.-F., Mackay 

D. , 2002, MNRAS, 336, 1364 

Jeffers S., Barnes J., Cameron A., Donati J.-F., 2006, MN- 
RAS, 366, 667 

Johns-Kruh C, 1996, A&A, 306, 803 

Johns-Kruh C, 2007, ApJ, in press (astro-ph/0704.2923) 

Johns-Krull C, Valenti J., Hatzes A., Kanaan A., 1999, 
ApJ, 510, L41 

Johns-Kruh C, Valenti J., Koresko C, 1999, ApJ, 516, 900 
Konigl A., 1991, ApJ, 370, L39 
Kiiker M., Riidiger G., 1997, A&A, 328, 253 
Kurosawa R., Harries T., Symington N., 2006, MNRAS, 
370, 580 

Kurucz R., 1993, CDROM # 13 (ATLAS9 atmospheric 
models) and # 18 (ATLAS9 and SYNTHE routines, spec- 
tral line database). Smithsonian Astrophysical Observa- 
tory, Washington D.C. 

Landi Degl'Innocenti E., 2004, Polarisation in Spectral 
Lines. Vol. 307 of Astrophysics and Space Science Library, 
Kluwer Academic Publishers 

Matt S., Pudritz R., 2004, ApJ, 607, L43 

Mohanty S., Jayawardhana R., Basri G., 2005, ApJ, 626, 
498 

Padgett D., 1996, ApJ, 471, 847 

Petit P., Donati J.-F., Auriere M., Landstreet J., Lignieres 

F. , Marsden S., Mouillet D., Paletou F., Toque N., Wade 

G. , 2006, in Casini R., Lites B., eds. Solar Polarization 
Workshop #4 Vol. 358 of ASP Conf. Proc. 

Petit P., Donati J.-F., Cameron A., 2002, MNRAS, 334, 
374 

Romanova M., Utsyugova G., Koldoba A., Lovelace R., 
2004, ApJ, 610, 920 



© 0000 RAS, MNRAS 000, 000-000 



18 J.-F. Donati et al 



Romanova M., Utsyugova G., Koldoba A., Wick J., 

Lovelace R., 2003, ApJ, 595, 1009 
Safier R, 1998, ApJ, 494, 336 
Schevchenko V., Herbst W., 1998, A J, 116, 1419 
Shu F., Najita J., Ostriker E., Wilkin F., Ruden S., Lizano 

S, 1994, ApJ, 429, 781 
Siess L., Dufour E., Forestini M., 2000, A&A, 358, 593 
Smith M., 1994, A&A, 287, 523 

Symington N., Harries T., Kurosawa R., 2005b, MNRAS, 
356, 1489 

Symington N., Harries T., Kurosawa R., Nay lor T., 2005a, 

MNRAS, 358, 977 
Valenti J., Johns-KruU C, 2004, Ap&SS, 292, 619 
von Rekowski B., Brandenburg A., 2004, A&A, 420, 17 
Yang H., Johns-Krull C, Valenti J., 2005, ApJ, 635, 466 



© 0000 RAS, MNRAS 000, 000-000 



